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Grain growth experiments on thin metalUc films have shown the geometric and topological char- 
acteristics of the grain structure to be universal and independent of many experimental conditions. 
The universal size distribution, however, is found to differ both qualitatively and quantitatively 
from the standard Mullins curvature driven model of grain growth; with the experiments exhibiting 
an excess of small grains (termed an "ear") and an excess of very large grains (termed a "tail") 
compared with the model. While a plethora of extensions of the Mullins model have been proposed 
to explain these characteristics, none have been successful. In this work, large scale simulations of a 
model that resolves the atomic scale on diffusive time scales, the phase field crystal model, is used 
to examine the complex phenomena of grain growth. The results are in remarkable agreement with 
the experimental results, recovering the characteristic "ear" and "tail" features of the experimen- 
tal grain size distribution. The simulations also indicate that while the geometric and topological 
characteristics are universal, the dynamic growth exponent is not. 
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Most metals, ceramics and minerals are polycrystalline 
materials containing grains of different crystal orienta- 
tion. The size, shapes and arrangments of these grains 
strongly affect macroscale material properties, such as 
fracture, yield stress, coercivity and conductivity. In 
magnetic systems, for example, the coercivity (or mag- 
netic 'hardness') can change by four or five orders of 
magnitude with a change in grain size [1]. Thus, un- 
derstanding and controlling polycrystalline structures is 
of great importance in the production of many engineer- 
ing materials and has motivated numerous experimental 
and theoretical studies of grain growth. 




FIG. 1: Schematic comparison between atomistic description 
of polycrystalline material and coarse grained picture of a 
smooth grain boundary network. Shown is a low angle grain 
boundary with aligned dislocations and two high angle grain 
boundaries in an otherwise hexagonal lattice. 

Grain growth in thin metallic films is one example 
where extensive research has been conducted. One very 
interesting experimental finding in such systems is that 
the grain size distributions and topological characteristics 
appear to be independent of many experimental condi- 
tions ^. More specifically, it has been found that for a 



large collection of Al and Cu thin films a universal grain 
size distribution emerges that is independent of the sub- 
strate, annealing temperature, purity, thickness and an- 
nealing time. Unfortunately the universal distribution 
is qualitatively and quantitatively different from the re- 
sults of extensive computational studies on grain growth 
[3l S] , which are based on the Mullins model [5 . In this 
model the problem is reduced to the evolution of a grain 
boundary network by relating the normal velocity Vn to 
the curvature n of the grain boundary, Vn = /i7A^, with 
mobility ja and surface tension 7, and specifying the Her- 
ring condition ^ at triple junctions. Various attempts 
have been made to extend the Mullins model and to in- 
clude more realistic effects, such as interactions of the 
film with the substrate, anisotropy in the grain bound- 
ary energy and mobility, grain boundary grooving, and 
solute and triple junction drag. However, no single cause 
is able to explain the experimental grain growth behavior 

mil]. 

In addition to the grain size distribution, the rate of 
growth of the average grain size has also been examined 
in detail. The Mullins model and its extensions all seem 
to predict that the average grain size, represented by its 
radius r{t) has a power law behavior of the form ~ t^/^, 
which follows immediately from the linear relationship 
between grain boundary velocity and curvature. Experi- 
mentally a much slower coarsening or even stagnation of 
grain growth in thin films is observed. This may be due 
to the fact that the Mullins model ignores the crystalline 
structure of the grains, the dissipation due to lattice de- 
formations and the Peierls barriers for dislocation mo- 
tion. It is difficult to reconcile Mullins type models with 
the atomistic features of grain boundaries, which (for low 
angles) can be seen as an alignment of dislocations where 




FIG. 2: a) - c) Grain structure obtained from postprocessing a PFC simulation at an intermediate time. The color coding 
indicates the averaged local lattice orientation for each of the maxima in the density field. An enlargement by a factor of four 
is used for each figure. The grain boundaries are visualized, d) Visualization of the density field, where each maximum is 
associated with a dot. The red lines indicate the grain boundaries. 



the driving force for grain growth is the stress associated 
with dislocation motion. The differences of the descrip- 
tion are shown schematically in Fig. [l] 

Atomistic descriptions can incorporate the important 
physical features missing in the Mullins model and have 
led to some important observations. It has been shown 
that the complex dislocation structure along curved grain 
boundaries gives rise to a misorientation-dependent mo- 
bility ^. Further studies indicate that grain boundaries 
undergo thermal roughening associated with an abrupt 
mobility change, leading to smooth (fast) and rough 
(slow) boundaries [9 , which can eventually lead to stag- 
nation of the growth process. The defect structure at 
triple junctions can lead to a sufficiently small mobility 
limiting the rate of grain boundary migration [TQI [TT] . 
Also tangential motion of the lattices are possible. For 
low-angle grain boundaries, normal and tangential mo- 
tion are strongly coupled as a result of the geometric 
constraint that the lattices of two crystals change con- 
tinuously across the interface while the grain boundary 
moves [T2j . As a consequence of this coupling, grains ro- 
tate as they shrink which leads to an increase in the grain 
boundary energy per unit length, although the overall en- 
ergy decreases since the size of the boundary decreases 
p!3HT5] . Each of these phenomena can be simulated using 
molecular dynamics (MD), see [16] for a review. How- 
ever, to study the effect of these phenomena on scaling 
laws, grain size distributions or stagnation of growth re- 
quires a method which operates on diffusive time scales. 

In this work, we employ the phase field crystal (PFC) 
method pT]. The model has been shown to successfully 
model grain boundary energies as a function of misorien- 
tation [18 and non-classical grain rotation during grain 
shrinkage and drag of triple junctions [T9]. In addi- 
tion lower coarsening exponents were already observed 
for hexagonal lattices [2QH22] . The aim of this Let- 
ter is to use the PFC model on large scales to obtain 
statistical data for scaling laws and grain size distribu- 
tions and to compare them with simulation results of the 
Mullins model and experimental data for thin metallic 



films. Since the experimental results in [2 seem to be 
universal, we do not fit the PFC parameters to a spe- 
cific material but consider an artificial setting within the 
simplest PFC model introduced in [17 . In dimensionless 
form the equation reads dtil^ — rV^|^, where the or- 
der parameter V^, is related to the time- aver aged atomic 
density, t is time, F is the mobility and free energy is 
T = /V^(-e + (V^ + l)¥/2 + V^V4) dr, see [H [23]- 
25 for details on the relationship to classical dynamic 
density functional theory (DDFT) and material specific 
parameterization. 

Fig. [2] shows a snapshot of a typical simulation. All 
simulations are performed in a periodic domain of square 
size L = 8, 192 starting from a randomly perturbed con- 
stant value of the particle density i\). After an initiation 
phase in which the white noise is damped rapidly, grains 
nucleate, grow and impinge on one another. Thereafter 
the number of maxima in the particle density '^ remains 
mainly constant and coarsening starts. Statistical results 
are collected after grains have reached a minimal size of 
100 atoms. Fig. [3] shows the obtained scaling results for 
the average domain area as a power law in times, i.e., 
t^, where q is 1/2 in the Mullins curvature driven model. 
In our simulations it is not clear that this relationship 
is valid as the value of q can be see to change in time 
and be dependent on the parameters of the simulation 
and initial conditions. For case "A", we either obtain 
an initial value oi q = 1/3, which turns into q = 1/5, 
or a constant value of q = 1/5, depending on the initial 
grain size. The constant scaling exponent is observed for 
larger initial grains. For case "B", corresponding to a 
softer material, the growth exponent increases to a value 
of q = 2/5, whereas for case "C", a harder material it 
decreases to q = 1/20. For all three cases, the growth 
exponent is significantly lower than the expected value 
(7 = 1/2 for the Mullins model. Similar low coarsen- 
ing exponents have been found for hexagonal lattices in 
[20] I2T] and in experiments for thin films of CoPt and 
FePt \2^. Extensive computational studies in [211 [22] 
further show a strong dependency of the scaling expo- 



nent on additional noise, which enhances the coarsening 
process. It has also been noted [22 that the addition 
of higher order time derivatives can change the growth 
exponent, which may be appropriate for three dimen- 
sional samples. In two-dimensional thin films (i.e., films 
with columnar grain structures), however, it is expected 
that the substrate/film coupling provides an effective fric- 
tion for rotation or translation that eliminates the need 
for such corrections. In either case, it is likely that the 
growth exponents are transient, because for very large 
gain sizes the Peierls-Nabarro barriers are likely to in- 
hibit further coarsening. This effect already occur at 
early times for quenches to lower temperatures, as con- 
firmed for points in the phase diagram in the solid region 
at (V^o,e) = (-0.31,-0.25) and (V^o,e) = (-0.29,-0.18) 
which show a frozen configuration. 
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FIG. 3: Mean area as a function of time together with the 
fitted scaling exponents for various points in the phase di- 
agram depicted in the inset. "Al" and "A2" have differ- 
ent initial grain sizes (Al < A2), the parameters are "A": 
(V;o,e) = (-0.29,-0.25); "B": (V^o,e) = (-0.25,-0.18); "C": 
(^o,e) = (-0.31,-0.30). 

While it is not entirely clear if there is a single, well- 
established dynamical exponent, the grain size distribu- 
tion functions appear to be much more robust. Fig. |4] 
shows the grain size distributions of the PFC simulations 
for the considered points in the phase diagram together 
with the experimental results of ^ and the results of the 
Mullins model taken from [3]. 

A considerable discrepancy between the experimental 
results and the Mullins model is already discussed in P, 
[7]. They differ in two important respects. First, the 
experimental grain structures have a larger number of 
small grains as evidenced by the peak of the experimental 
reduced area probability density residing to the left of 
that for the simulations based on the Mullins model, a 
feature that has been termed the "ear". Second, the 
experimental grain structures have "tails" that extend 
to significantly larger sizes than those seen in simulations 
based on the Mullins model. While only very few grains 




FIG. 4: Grain size distribution with reference to radius (area 
in inlet). Shown is the mean distribution, obtained as the 
average of the last time steps in the considered cases in Fig. 
|3] The curve are fitted to a lognormal distribution with pa- 
rameters (/J^jCr) = (—0.13,0.53). The experimental data and 
the results of the Mullins model are taken from 2 



seen in simulations exceed 4 times (and only rarely do 
they exceed 5 times) the mean area, the experimental 
grain structures exhibit maximum grain areas that are 
between 8 and 42 times the mean, with a sizable fraction 
of grains whose areas exceed 4 times the mean grain area 
(~ 3% by number, representing ~ 18% of the total area). 

Various closed form distributions have been proposed 
to fit the results of the Mullins model, e.g. the Louat, 
Hillert, Rios and Weibull distribution (see [3 and the 
references therein). They all not only differ in the "ear" 
and "tail" region, but they also peak at r/{r) > 1, again 
in disagreement with the experimental results. The PFC 
simulations not only recover the qualitative behaviour 
of the experimental results, they almost perfectly fit the 
distribution, and can be very well described by a log- 
normal distribution. The grain size distribution appears 
to be self-similar which is analysed in detail for case 
"Al" in Fig. [5] All results are obtained without ad- 
ditional noise. However, simulations that included noise 
(not shown) produced grain distributions consistent with 
the zero noise case. Further analysis indicates that, also 
in agreement with the experimental data, small grains 
are primarily 3 and 4-sided, whereas large grains have 
primarily more than 6 sides. 

The importance and prevalence of the formation and 
properties of polycrystalline materials has lead to an 
enormous amount of theoretical and experimental re- 
search. Unfortunately theoretical progress has been hin- 
dered by the lack of computational methods that can cap- 
ture the essential physics on the time and lengths that 
are appropriate for such phenomena. While MD simula- 
tions are currently unable to reach time scales required 
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FIG. 5: a) Grain size distribution with reference to radius 
at the labled times in the inlet, cooresponding to case "Al" 
in Fig. [3] in comparison with the experimental results from 
[2] . The initially narrow distribution broadens rapidly and its 
peak shifts towards smaller grains. For large times the grain 
size distribution appears to be self-similar which is further 
illustrated in b) showing the time evolution of the parameters 
G and /i of a log-normal distribution fitted to the considered 
snapshots, again in comparison with the experimental results 
from J2J shown as the horizontal solid lines. 



to observe self-similar growth regimes, coarse grained de- 
scriptions based on the Mullins model seem to lack the 
essential atomistic features allowing for bulk dissipation 
during grain growth. In this work large scale numeri- 
cal simulations of the PFC model were used to exam- 
ine the phenomenon of grain growth in two dimensional 
systems. The results of these simulations are in remark- 
able agreement with universal aspects of the geometric 
and topological characteristics of the grain structures in 
thin metallic films. Among other features they capture 
both the "ear" and "tail" characteristics of grain distri- 
butions that have proven difficult to obtain with previous 
models and methods. Thus the PFC model provides a 
key resource for future research in which realistic grain 
structures are required. Although not examined in this 
work, the model also incorporates mechanical properties 
of the system and thus can be used to study, for exam- 



ple, the relationship between growth conditions and the 
structural stability of polycrystalline materials. 
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